clear

[point, seg, tri]=importfilemesh('Th.msh');
u=importfiledata('phi0.bb');
% ux=importfiledata('ux.bb');
% uy=importfiledata('uy.bb');

% pdeplot(point,seg,tri);
pdeplot(point,seg,tri,XYData=u, ...
        mesh="on", ...
        colormap=[0 0 1; 1 0 0]);
% pdeplot(point,seg,tri,XYData=uy,...
%         mesh="on");



